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ABSTRACT 


A luathejnatical roprosentation has been clcveXopecl fob the 
electromagnetic force field, the fluid flow field, the 
temperature field (and for transport controlled kinetics) , in a 
levitation molted metal droplet. The technique of mutual 
inductances was employed for the calculation of the electro- 
magnetic force field, while the tur):>ulcnt Navicr - Stokes 
, equations and the turbulent convective transport equations 

were used to represent the fluid flow field, the temperature 
% 

field and the concentration field. The governing differential 
equations, written in spherical coordinates, were solved numerically. 
The computed results were found to be in good agreement with 
moasurementj reported in the literature, ’‘"yarding the lifting 
force, and the average temperature of the specimen and carburi?ia- 
tion rates, which were transport controlled. 


1. INTRODUCTION 


The levitation melting technique, that is the positioning of 
metal droplets by an electromagnetic force field, generated 
by induction coils of suitable geometry ha^ gained a great deal 
of popularity during the 1950 's and 1960’s. The principal 
attraction of this technique v;as that small metallic samples (say 
ranging in weight from .5 to 5 g.) could be melted and reacted in the 
absence of a solid container. It was generally held that the eddy 
currents induced in the metallic sample produced 'a spatially uniform 
temperature field thus the technique was foynd to be particularly 
useful for studying gas - melt interactions. 

More recently levitation melting is gaining rather greater 
prominence, because electromagnetic force fields are being used 
for the positioning of samples in space processing experiments. 

One of the important features of these systems is that the . 

absence of gravity will necessarily eliminate thermal 
buoyancy driven natural convection, a factor, of particular 
importance in certain solidification studies. However, since the 
samples are positioned by an electromagnetic force field, the question 
will naturally arise whether this force field in itself may give 
rise to convection currents, furthermore whether any temperature 
gradients that may exist on the free surface of these samples 
may give rise surface tension driven flov;s. 


ClcarJy, a quanUitaUivo undorstanclinq of Uiosc phenomena would 

be highly desirable, for both the planning and the interpretation 

of in-flight experiments. Furthermore, the quantitative • 

description of the temperature and velocity fields in levitated 

droplets under conventional ground based experimental conditions 

would also be of considerable interest, since this knowledge 

would be very helpful in the interpretation of any kinetic measure- 

» 

ments. 

Previous work in the area of levitation molting may be 
divided into two categories, namely experimental studies concerned 
with physico-chemical measurements at high temperatures, which 
did not consider transport phenomena in the melt on the one hand, 
and purely elcctrodynamic work, concerned with the calculation of 
the lifting force on the other. 

The principal purpose of the present paper is to develop the 
theoretical framework for combining the elcctrodynamic and the 
transport aspects of this oroblem. 

2. ITORMULATION 

Fig. 1 shows a sketch of the conventionally used coil arrangonents 
for the levitation melting of metal droplets. It is seen that an AC 
current is being passed through two (or more) opposing coils and as 
a result an electromagnetic force field is being generated. 

This force field will have two effects on a metallic specimen • 
located between the coils: 


(i) A net foirce (lifting fojrco) will be exerted on the 
specimen, which, when large enough compared to the weight of 

the sample, will be able to balance the gravitational force and thus 
"levitate the specimen". 

(ii) The induced eddy currents dissipated in the sample will cause 
"Joule Heating" and will also generate a stirring force. 

In developing a quantitative representation of such a ] 

* * * 

system, one has to 

(a) Calculate the electromagnetic force field acting on the 
specimen, hence the lifting force, i.e. the conditions for levitation 
for a given applied magnetic field, such as the coil current and the 
coil confiauration. 

(b) Calculate the fluid flow field, as generated by the 
electromagnetic (Lorentz) forces. 

(c) Calculate the temperature field that x'esults from the 
balance between heat generation and heat dissipation due to 
convection and radiation. 

Here one should remark that tlie governing equations are necessarily 
coupled, because the fluid flow within the sphere is driven by 
the combined effects of thermal natural convection and electro- 
magnetic forces. Furthermore, the spatially variable heat 
generation in the system is caused by the current distribution. 

While a number of electromagnetically driven flow problems have 
been successfully modelled in the past,' , the present system poses 
the additions complication of spherical^ symmetry. 


£,1 _Calculation^of Uie KlectronmgnGHc Potco Field ancl 

I 

of the Lifting Force 

Let us consider a molten sphere having spatially uniform 
» 

physical properties sketched in Fig. 2. Lot us consider 
furthermore, that the applied electromagnetic force field also 
shown in the figure will be axially symmetrical. For spherical 
geometry the induced current will be flowing in the (j) direction. 

The induced current will generate a magnetic field and the 
interaction of this magnetic field with the current will give 
rise to an electromagnetic force field. 

These electric phenomena will have three manifestations: 

(a) A lifting force will be generated 

(b) Thermal energy will be dissipated in the 
metallic speciinan 

(c) A stirring force will be exerted on the s?Decimen. 

F„ , the total levitation force is given by 

i 

; ' F. = / (1/2) Re (JxBp dv (1)** 

/ . ' 
t 

furthermore the power adsorbed by the melt (i.e. the total 
Joule Meat generated) is: . 

P = /(l/2)Re(J . J*)/o dv ^2) 

V 

finally the spatial distribution of the time-averaged stirring force in 
the melt per unit volume is: 

F = 1/2) Re (J X B *) 

S ' ♦ ^ 

'■*Thc list of symbols is given at the end of the paper. 


( 3 ) 


I 
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What necrls to be done here is to evaluate the induced current 
density J and the magnetic ilux intensity B, for a given coil 
geometry and coil current. 

in essence there are two ways in which this may be accomnlishod . 
One of those would involve starting with the differential Maxwell's 
equations, solving for the magnetic field and then computing 
the induced current. . 

The alternative technique, which will be .used here, involves 

I 

the concept of mutue\l inductances. In this integral formulation 
the soecimen (i.e. the levitated droplet) is considered to consist 
of a set of electrical circuits and one. solves for 
the current generated in each of these. Once the current 
distribution is known, the magnetic flux density is readily calculated 
and hence the terms appearing in Eqs. (l)-(3) may be evaluated. 

The magnetic flux density B for a given current distribution 
is given by Arapere's law: 


dxa. 


B = -j 


‘4 13 


/ - 
vol r 


,2 


dv 


(4) 


where J is the current density in a volume element dv at a distance 
r' from the point at which the flux density B is being evaluated 
and a is the unit vector in the r direction. The quantity J 

17 *si» 

must include the induced current within the levitated sphere as 
well as the current within the coil. 

Let us define the vector potential A, as; 

t ‘ ' 

B = V X A • (5) 


Then Eq. (4) becomes 

J 


o 

4 II 




dv 


( 6 ) 


vol r 

On noting the axial symmetry of the field, the levitated 


H. ' • V 

Sphere is divided into elementary circuits as shov;n 

in Fig. (3).‘‘Rach segment of the spherical shell is considered 



as a region of constant current density. For such a system 
Eq. (6) may be written as; * . 



^^0 

I I 1 11 . 1 — 


sphere • dJ-„ coil^ * 


(7) 


where dJl is the line element of circuit of constant current density, 

, » ' 

and S is the' cross sectional area of the circuity and I(k) is the coil 
current 


It should be noted that first term on the rhs of Eq. (7) 
describes the induced potential while the second term describes the 

applied potential. 

Prom Faraday's law of induction, the induced current in anv 
circuit, i in a levitated sphere is given explicitly by; 


E. .dZ. s. - ^ / B.ds 
1 X dt „ 


( 8 ) 


where the surface of integration is bounded by the periphery of the circuit. 
For’, a time harmonic field we have 


J = Re { J 
B = Re { .B 


(9) 

( 10 ) 


Using Eq. (9) - (10), Eq. (8) can be written as: 

SDhere n 

0 J. . df. =OwO { T. M. + I M. I(k)) 


( 11 ) 


1, c 


^X *^'C 


where 


u 

^00 

4 IT 


( 12 ) 


4 
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io tho nmfciml inductance, (on seli inductance, when i»C) which may 
be calculated using the ring-iring formulas (8). 

Finally in order to avoid dealing with complex numbers, equation 
(12) is separated into real and imaginary parts, as follows: 

The levitated sphere depicted in Pig. (B) is divided into 
N X L circuits. The eddy current distribution can be obtained by 
solving 2 x N x L real, simultaneous algebraic equations. 


ll * * 

coils 

Ll "(i,« k- 


N L 

j: 

vi=l V=1 




k=l 


m 


In the present case, the vector potential has only one 
component, in the <|) direction; since the vector potential has the 
same direction as the current, B^, and are x'elatod to A as follows 


B “ TTT (A.iSi.n 0) 

,r r sin u c't) (j) 


(15) 


B. 


r Ir 


^ <lf') 

Thus finally the knowledge of B^, B^ and enables one to evaluate 
the quantities P^ , P and which were defined in Eqs. (1-3). 


2.2 Calculation of tho Fluid Flow Field 

In gcnoral the fluid flow field in levitated metal djropletn will 
be driven by two forces, that is the electromagnetic force field 
(which we have just calculated) and the buoyanev force field, which 
is due to temperature differences v;ithin tho specimen. It may bo 

I 

expected that tho velocity field will be turbulent, so that aonrooriato 
techniques have to be invoked for evaluating the turbulent shear 
stresses (i.e. the Reynolds stresses). 

The following principal assumptions arc maclo* in the development 
of the governing equations. 

(1) Tho flov; field will be considered to be two dimensional, 
that is only the r and the G components of the velocity will appear. 

(2) The k-F. model^^'^^^ may be used to represent the turbulent viscosity 

(3) Anv possible damping of the turbulence by the magnetic field 
is neglected 

(4) No allowance is made either for the anisotroj.-^ of tiir- 
bulence near the free surface, or for the possible damoing effect due 
to surface tension forces near the free surface. 

On commenting on the appropriateness of these assumptions it is 
noted that because of the axial symmetry (symmetrical location of 
the specimen regarding the coils) the assumption concernincr the two 
dimensional nature of the velocity field was thought to be reasonable. 

The k-e model has been selected because this has been a quite 

widely used approach for modelling turbulent electromagnetically 
(11 12 ) 

driven flows , v/ith considerable success and at present 

there appears to be no obvious, more satisfactory alternative. 


‘ K 
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'jPhc Iftst Uwo assumptions repirosont an ovnx*Bimplif ica- 
tion# which ouqhfe to be refined, once experimental data become 
available in order to permit this. The Implications of these 
assumptions, in the light of some experimental measurements, will 
fee discussed subsequently. 

For steady two dimensional flow in spherical coordinates, the 
governing equations take the following form; 

The equation of continuitv; 


J2 h r-sIHT Iff “ » 

The equation of motion j 
r direction 

1 


(17) 




^f£ + 


1 ^ 

^ (py u sin 0 
rr r srn 0 90 ’-'r 0 


- n -- sin 0 S rt\ cri 
eff rO) ~ SUj. 


(18) 


0 direction 


ij Ip <P!Jj.Uo r'siTnr” Iff *'’'"3 ° 


^‘eff ^ ^^0 


U9) 


where S. . is the symmetric part of the strain rate tensor 
r 3 


SU^ - 


2 

’JO 




- If - Vf 0 


+ i pRe ( B*) 


( 20 ) 
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^ X ^0 1 r>»u 

SUg 


of£^4.(|. 


+ '(n-p^^) geos 0 + I ^ \ 


whore 


P * P *" P^ gi^ <208 0 

p » p^ (i - a v^r) 


®rr “ 


S * i. . 0 + ^ 

*=*00 r (TIP r 


c- « U. UgCOt 0 
^A{|, ® -£ 4- -il — 

t ■ r 


1 a 1 

■I lr< ? -^> 


The turbulent or eddy viscosity may be calculated 
by using the k-c model: 


P kVe 


Separate transport equations will have to be written for k, the 

turbulent kinetic energy and c the turbulent kinetic energy 

dissipation; thesg^talce the following form; 

* * * » « 

* 

The turbulent kinetic energy; 


r 


2, ^ ^'t , 9kv 

+ o7> 


1 9 

r sin 6 W 




(PU^ sin Ok- (h t 7j|-) 


sin G 9k 
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where 

®\ - 2»l: “"L + ®00 + ®r0> - 

The turbulent energy dissipation; 


(30) 


1 

•~2' 

r 


8 

•gr 


(Ptj^r^r. 




0G 


J)r.v 


d 

W 




(PUf, sin 0 e - (wH- 


sin 0 


l£- 

90" 


su. 


(33) 


Whore 


SU *s 2C, -r- 
1 k 


«00 K + ®r0> 


- Cj P 


k 


(32) 


'and Oj^f and are constants. At high 

• * I 

values of the Reynolds number all these quantities are generally 

aesumed to be constant. At low Reynolds numbers and C 2 will 

? * 

become dependent on Ue^K^c'/vs). The values of these constants and of 

(9,10) 

the func iqn are given cljsev’hore. ; 


2.3 Temperature field 

The temperature field in a levitated metal droplet is governed 
by the rate of heat generation, due to Joule Heating, the rate of 
convective and radiative loss from the free surface and by the rate 
at which thermal energy is being transported inside the specimen, 

4 

by conduction and ,by convection. 

Written in spherical coordinates the convective energy balance 
equation takes the following form; 


iwiiimii lit ur TTiiTi^' "• , 




~13~ 




3 * L.(n JL2. 
r 




U . '‘t . DT» 1 




^ sXnU 




Re fJ».7*) 


_______ <33) 

‘pr« Pr|,' 2 PC 

^ U Q 

The boundary conditions necessary to complete the statement of the 

* *• * * ♦.* 

nroblom arc cjivon as tollows: • •• * ' 

At the axis of symmetry 0 = 0 and 0 « ir 


U 


0 . 


0 ) 


3Uj730 = 0 
3K/30 =- 0 
3e/30 “ 0 
3T/90 - 0 


Y? ntatement cf symmetry. 




At the surface of sphere r = R 

3U^/Sr = 0 

^'r = 0 . 

3k/3r = 0 
3e/3r = 0 

- k = h + Of.(T -T^) 


(34) 


(35) 


These surface boundary conditions express the physical 
constraints that momentum, turbulent kinetic energy and the rate of 
turbulent energy dissipation are not transported across the phase 
boundary , and that the velocity normal to the phase boundary is zero 
The last equation expresses the continuity of the heat flux. 


h, the convective heat transfer coefficient between the levitated 

drop and the surroundihq qas stream was correlated for forced. 

f 13) 

convection and natural convection as ' 


Nu 0 . 81^0 Pr^/^ 

150< Rg< 500 

3d) 

Nu « 0.78 (Gr 

10^ < Gr < 10® 

■ (37)' 


3. COMPUTATIONAL TNC UNIQUE 

In generating a solution of the governing equations the 
magnetic field equations were solved first. This involved 
divd.ding the sphere into 8x15 circuit elements and solving the 
resultant set of linear equations. Once the current distribution 
was known, the magnetic flux density was readily calculated, 
together with the electromagnetic force field vector, the pattern 
of Joule Heat generation and the lifting force. 

In the subsequent phase of the computation the fluid flow 

equations, the convective heat flow equation and the conservation 

equations for k and e v/ere solved together, using a finite difference 

d'4")’ 

technique described by Spalding. ' 

A 10x17 grid structure v;as used and the computation required 
about 100 seconds on MIT's IBM 370 digital computer. 

4. ■ CQMPUTKD RESULT S 

In the following we shall present a selection of the computed 
results, pertaining to the lift force exerted on levitated droplets, 
on the fluid flow and heat flow phenomena and finallv this informa- 
tion will be used for the interpretation of some measurements, 

obtained by one of the authors on the carburization of molten iron 
(13,15) 


droplets. 
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The actual calculations v;ere carried out for input parameters 
which correspond to the experimental conditions used for these 
carburization measurements. One gram of iron was levitated in a 
CO 2 +CO gas mixture, at 40 atmospheres, with a gas flow rate of 
5 1/m (STP) . The tube containing the droplet and the gas stream 
was 13 mm in diameter . Under steady state conditions the temperature 
of the molten drop was 1650^0. 

The actual coil configuration used is given in Fig. 1; the 
coils were made of silver tubing, 3.1 mm in diameter. Two co- 
planar turns, 17.2 mm inner diameter were employed, both above 
and below the levitated sphere. The distance between the upper 
and the lower turns was 10 mm. 

• The pov/er supply was a lOkW, 450 kHz high frequency generator. 
The actual coil current v;as not measured, but its value will be 
deduced from subsequent heat flow calculations. 

4.1 Lift Force Calculations 

In order to find the position of the levitated sphere, for an 

* 

estimated coil current, the lifting force and pov/er absorbed is 

calculated along the axis of the coil. Figs. (4) to (6) show the 

results for coil currents of 200, 250, 300 Amp respectively. The 

symmetry of the force field and power absorption around the center 

of the coil is a result’ of the symmetry of coil. ^t a 200 Amp coil 

-3 

current. Fig. (4) , the maximum liftxng force is 7x10 N which 

-3 

is less than weight of the drop (9.81x10 n) . Under these 
conditions it is not possible to levitate the 1 gram iron drops. At 
higher coil currents Figs. (5) and (6) levitation is possible as 
indicated by the intersection of the curves representing the lift 
force (Fj) and the lA^ight of the particle (mg). There are two 
positions where the electromagnetic force balances the 


wolqht or the. (Iron, an cloarlv shown in PXq. (5) » but on3y Xho 
larcfor of t liose two vaUu>n dotni oorrcsnoml to a nlabXe equilibrium, 
sincG any upwaixl or clownwanl tl^viation from it produ('’cs a 
restoring lifting force tloorcment or increment reenoctivoi y . 

It is inl'orosting to note that at Koro gravity, the stable position 
corresponds to the center of the coil. The power absorbed at ^!50 
was 70 VJ while for 300 this figure was The decrease of the 

power x’^iliisorption by increasing coil current is the basic feature 

« 

of electrcrnairnetic levitation. • . 

The coil current may now be e.stimaled wilh reasonable accuracy, 
because a value of about 250A will give tlu' un*oper lifting force. 
Furthermore, as will be shown subsecmentlv, the corresnonding 
power innut will give the correct droplet temnerature, usin<) nj’< 


overall fsnergy balance. 

4“, 2 F. 1 e c t rom’a ii n e tic 


forc es a nd hefft 5ienprvU 


ion in IcA/^i tilled 


sohere. 

»CT :*a2-» 1 

Fi.cf. (7) shows the computed force field and the heat generation patten 
in the. drop. It is seen that the oleetromagnotic field is markedly 


attenuated 

I 


the skin denth which is 


0.9 mm for molten iron at 


450 KUz . As a result of positioning the levi.tatod snherc dose to 
to the bottom coils the force field and lieat generation are moved 


toward the lower half of the snherc. 


clron 

* 


.gince the current only penetrates a very small region in the 
the magnetic flux density inei*eases rapidly with r as it 


apnroaches to the surface, therefore the electromagnetic force,s are 
almost normal :^o‘ t!ie r direction. 
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^^3 F lov^? FA olds 

The computed velocity field and tempeirature field are given 

in Fig. (8) . The computed velocity field shows n double loop 

circulation pattern, which is consistent with the force field shown 

in the previous figure and is to be expected for a stationary 

* 

electromagnetic field. The maximum velocity at the free surface 
was found to be about 0.28 m/s, which is in reasonable, iqualitative 
agreement with measurements reported in the literature. More, 
specifically Robertson estimated surface velocities ranging from 
about 0.1 - 0,.2 m/s for a similar system considered here.^^.^V 

An important point to be made about the temperature field is that 
the average temperature of 1650^C agrees well with the measurements. 
The examination of the temperature field is instructive, 
because it shows that a temperature difference of about 10°C may 
exist even in this relative.ly small scale system, in spite of 
regions having quite high linear velocities. This is a significant 
finding which deserves further comment. 

The apparently quite rapid circulation rates that were predicted 
do not guarantee that the system is \ ull mixed, provided the time 
scale for convection is much smaller than the diffusive time scale. 

For the present case, t^, the convective time scale is given 
as ; 

t £! 0.02 s 

c U 

while the diffusive time scale (for turbulent flow) is given as: 
a 2 
“eff 


( 33 ) 




I 


.This can explain the existence ol; the slight temperature gradients 
within the system. The question of mixing is of greater imoortance 
in diffusion problems (since the molecular diffusivity for molten 
metals is much smaller than the thermal diffusivity, even v/hen 
both are. effectively enhanced by turbulence) and this question 
will be discussed subsequently. 

I 

Fig. 9 shows a computed plot of both the turbulent kinetic 

» 

energy and of the ratio; effective viscosity /iftolepular viscosity 
for the system. The important point to note here is that k is 
sharply reduced and that the ratio tends to 

a value of 2 on approaching the f,ree surface. It has 
to be stressed that this behavior has not been externally imposed 
on the system, but ratlier is a direct consequence of the curvi-linear 
geometry employed here. More specifically the rate of energy 
transfer from the center toward the outside is being diminished 
by the generation of angular momentum. Another important factor, 
which contributes to the laminarization of the flow field close 
to the free surface, is the energy consummed in the acceleration 
of the fluid in this region. " 

It should be noted, th^t yet another factor, which may 
contribute to the laminarxzation of the flow field in the vicinity 
of the free surf ace, is the role played by surface tension forces. 

* This effect has not been considered explicitly in this formulation, 

■ but will be commented upon in the subsequent section of this paper. 

* . t 

% 

n; 




% 

I 

imi'ifirii — ti-er n — 
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4.4 Reaction Kinot.lcs in a Levitated Dronlcfc 
'jCho previously given computed results were found to be in 
good quantitative or at least somi-cqualitative agreement with 
measurements, regarding the lift force, the surface velocities 
and the mean temperature. Examination of the computed temperature 
field has indicated that the system v;as not very v/ell mixed, in 

spite of the quite rapid circulation rates, but no direct experi- 

♦ 

mental evidence was available to support tliese predictions. 

Much better insight reoarding the extent of mixing in tlie system 

mav be obtained bv the re-examinntion of some kinetic measur'^ments 

(15) 

reported by El-Kaddah and Robertson. 

In this work the rate of carburisation of iron droplets was 
studied with a CO^/CO mixture, for experimental conditions 
which have been given previously. 

The carburi^.ation reaction involved the following reaction 
scheme. 

2C0 ^ CO 2 + (C) 

Previous studies have shown that the chemical reaction rate 
\ms very fast under the experimental conditions employed, further- 
more, the calculated gas phase mass transfer rate was found to be 
faster than the experimentally observed rates. It follows that 
the rate limiting step had to be liquid phase diffusion or a 
mixed control involving both the gaseous and the motion phases. 

On postulating mixed control the governing equations take 


the following form: 


-20 


coiiGorvation of the transferred species: 


i^(L- fin * O- r^* *<• -iL.) HJ£) 4. ^ 

yt r^'iir ^ 'Sc. Sc,/ ()r^ ^ r sTn 0 

/s V 




nc, 


1 

.szm^ 


3 /,,n n n sin 0 / n . 'U; » ?Cv ^ 

(pU^ SXn 0 C - H- rrjT) « 0 


0 r 'SCjj^ ■ Sc^/ 00' 

the boundary conditions are given as foiJ.ows: 


(4U 


C 0 at t “ 0 

specifying soro initial carbon content, and^ 


/ J 
'Sc" 


"t. 


■) 


ia 


=■• >^!I (>frn -X„o ) 

*y >jt>» fy 

B 


Sc^' 'r=R '''CO^-^CO 


(42) 


(43) 


specifying the continuity of the flux of carbon of the gas-melt 

interface where X ^ is calculated from equilibrium considerations and 

''•2 

s 

may be expressed as 


X 


CO, 


1 


K ^‘c 
’” 2 P 


\l 


Kfif* 

(It «2ir)'^ - 1 


where the activity of carbon is given by \ 
I’c 

= C 10 
c 


(18) 


(44) 


(45) 


This system of equations is readily solved, utilising the 
previously obtained values of the velocity field and of the 
turbulence parameters, needed to compute the values of the eddy 
dif fusivity . 

In view of the sharp changes in both the eddy diffusivity 

and in the volume of the grids in the vicinity of the free surface a 
somewhat finer grid structure w^as used in these calculations, 

employing 17x16 grids, with the extra po’ihus being located in the 
vicinity of the free surface. 
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Fig. 10 shows the experimental measurements with the full 
circles; the computed results # using turbulent flow considerations 
are given by the broken line, while the solid line shows the 
results of calculations, assuming laminar flow in the levitated 
sphere. 

It is seen that the assumption of laminar flow gave 

excollent agreement between the measuiements and the oredictions. 

The calculated results based on turbulent flow jnopcar to be 

> 

quite close also, but this may be somewhat misleading, because 
the curve drawn with the broken lino closely approaches the limit 
of gas phase mass transfer control. 

These results clearly indicate substantial damping of the 
turbulence in the vicinity of the free surface, v;hich appears 
to go well beyond the behavior that was attributable to the 
acceleration of the flow and the curvi- linear coordinate system 

I 

that were discussed earlier. 

Mechanisms that may be responsible for the damping of 

turbulence in the vicinity of the free surface may include the role 

played by surface tension forces, as discussed by LciVich^^^^ and Davis 

However another important mechanism that may have to be considered 

is the damping of the turbulence by the electromagnetic field, 

which would be particularly important in damping the velocity 

( 21 ) 

components parallel to the free surface. . . 

In view of the quasi-laminar behavior postulated in the 
vicinity of the free surface it may be of interest to show one 
additional figixre, depicting the computed velocity field and 
temperature field for laminar flow conditions. This is done in 
Fig. 11 where it is seen that the behavior is qualitatively similar 
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t.o that founcl foj: turbulent flow, except for the fact that the 
linciir velocities are some 30-40% higher and the isotherms arc 
consistent with a rather well defined laminar circulation pattern. 

5. CONCL UDI NG REMARKS 

A mathematical representation has boon developed for the 
electromagnetic force field, the fluid flow field, temperature 
field and transport controlled kinetics in a levitated metal 
droplet. 

, I 

The statement of the problem included the use of mutual 
inductances to calculate the. current density distribution in 
the droplet resulting from a given coil currentand coil configura- 
tion. Knowledge of the current density distribution enabled the 
calculation of the lifting force, of the stirring force within 
the sphere and of the heat generation pattern. 

The fluid flow field in the sphere was modelled by writing 

the turbulent Navior-Stokes equations , in spherical coordinates, with . 
symmetry about the angle 0. The effect of both buoyancy and 
electromagnetic’ forces was included in the formulation,,, thus the 
fluid flow equations v/ere coupled with the differential thermal 
energy balance equation, which included convection, eddy conduction 
and heat generation. A differential component balance was also 
developed to represent the carburization of an iron sphere, due 
to the decomposition of CO at the metal surface. The governing equations 
were solved numerically, using the k-e model for the turbulent 
viscosity and employing primitive variables. 

Computed results were obtained for the lifting force, fer the 
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heat gonorntion pattoimr tho velocity and temporaturo iJields 
and the rate of carburi^sation. Tho theoretical predictions 
were compared with experimental measurements., reported previously. 

In a macroscopic sense tho theoretically predicted lifting 
force and moan bulk temperature of the levitated sphere were 
found to be in good agreement with the measuromonts. 

It is thought, however that the detailed insight into the 
microscopic behavior of tho system, provided by these calcula- 
tions was perhaps of greater significance. 

The calculations have shown that the electromagnetic force 
field in the molten metal sphere produced quite high melt 
velocities, of the order of 0.1 - 0.3 m/s; these values were 
in qualitative agreement with the estimates deduced from the 
visual observation of levitated specimens. V7hilc the fluid 
flow field was definitely turbulent in the central part of the 
sphere, notwithstanding the quite small linear scale of the system, 
laminarization had occurred to a considerable extent in the vicinity 
of the free surface. 

In a mathematical sense this laminarization is attributable 
to the acceleration of the fluid An a curvi-linear system. In a 
physical sense laminarization near the free surface may also have 
been promoted by surface tension forces. 

This laminarization of the flow field has important consequences 
regarding other transport phenomena. The calculated temperature 

t 

fields within the sphere have shov;n that a certain non-uniformity 
of the temperatures may exist within the system (say temperature 
differences of up to 10°c in the particular case calculated) 
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notwithstnncling the high linear fluid volociuios and the relatively 
small physical size of the specimen ( 6mm diameter) . 

This lamiharization has an oven greater influence on transport 
controlled kinetics. In this regard the theoretical predictions 
for the rate of carburization were comoared with experimental 
measurements. The theoretical prodictionB were based on the 
assumption that gas phase mass transfer and liquid phase convective 
Uncluding turbulent) diffusion being the rate' controlling factors. 

It was found that the use of a turbulent model would have over- 
predicted tlie carburization rates, while the postulate of laminar 
flow gave very good agreement between mensurements and predictions. 

This latter finding clearly confirmed the important laminarizing 
effect of surface tension in the vicinity of the free surface, which 
in the present case corresponded to the major part of the physical 
domain. 

At this stage it ought to be mentioned that two assumptions 

implicitly made in the formulation, namely spherical geometry and the 

absence of surface tension driven flows may be justified a posteriori 

by the reasonable agreement between the measurements and the predicted 

behavior. The close approximation to spherical geometry has been 

observed by many investigators, furthermore, the absence of significant 

temperature gradients on the surface should preclude significant 

surface tension driven flows, especially as the electromagnetic forces 

are at maximum at the free surface. 

$ * 

In conclusion it may be stated that a, realistic mathematical 
representation has been developed for the electromagnetic force field 
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boat, mass and momentum transfer in levitated metal droolets and 
that the theoretical predictions appear to bo in very good agreement 
with measuromonts . 

It is suggested that these results may gain useful application 

I 

in two distinct areas. One of these is the rational planning and 

t 

interpretation of space processing cHperimonts, which involve 

levitation melt.ing. The other orovidos a potentially very useful 

» 

tool for the interpretation of gicund based levitation experiments, 
through the quantitative definition of the fluid flow field in 
these systems. 

Prom the standpoint of basic fluid dynamics the laminarization 
of the flow field in the vicinitv of the free surface is a very 
interesting finding, which will deserve further detailed study. 
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hIS'P OF SYMBO LS 

• MMWnnNMnMiniMMMMK Jtt- 

a^. activity of carbon 
o 

A vector potential (Wb/m) 

2 

B magnetic flux density (Wb/m ) 

m 

2 

B^ applied magnetic flux density (Wb/m ) 
c carbon concentration (Wtl) 

mm 

3 

C concentration diffusing species (kg/m ) 

C specific boat of levitated sphere (J/kg*^K) 

C^fC^fC^ constants of k-e model (m) 

t 

^5 electric field (v/m) 

«w. 

F electromagnetic lifting force* (N) 

% « 

2 

g gvMvitional acceleration (m/s ) 

3 

F„ eiectromagnetic stirring force (N/m ) 

•- b , 

2 o 

h heat transfer coefficient WW/m K) 

I coil current (A) 

Im imaginary value of 

2 

J induced eddy current density (A/m ) 

~ 2 2 
k turbulent kinetrc energy (ms) 

kg .n,ass transfer coefncient ■ . . 

in mass of levitated drop (kg) 

M. „ mutual inductance (II) levitated drop (W) 
p power adsorbed by Eg. 22) 

P pressure (N/m ) 

P modified pressure 

$ 

r radial coordinate (m) 

Be real value of 

# 

2 

S cross sectional area of elementary circuit (m ) 

T temperature (°K) 

Tg ambiant gas temperature (^K) 
t time (s) 

y velocity (m/s) 


I 

\ 

j 

i 


I 



t 


4 

X moXe fraction of reactive gases 

magnetic permeability (U/m) 

turbulent viscosity (kg/m.s) 

. 3 

p density of levitated sphere (kg/m ) 

2 3 

e turbulent kinetic energy (ms) 

C emissivity of levitated sphere 

a electric conductivity of levitated droplet (J2m) 

■ 2 o 4 

0 . Stefan-Boltzman constant (J/s.m K ) 

0, 0 constants in k-£ model • . 

k s ^ 

V kinematic viscosity (m'/s) 

(0 frequency 

^eff effective thermal diffusivity 

a thermal expansion coefficient 

Dimensionless group 

Nu Nussett number 

Pr Prandtl number 

Re Reynolds number 

Re,p turbulent Reynolds number 

Sc Schmidt number 


Gr 


Grashof number 


Figure Captions 

1. Sketch of the levitation coil used in this work. 

2. Sketch of the electromagnetic field in levitated sphere, 

3. Illustration of the elementary circuit used in 
numerical computation. 

4. Computed lifting force and power adsorption for a 

molten iron droplet along the axis of the coil. Coil 
current. 200 Amp. ’ 

5. Computed lifting force and power adsorption., for a 
molten iron droplet along the axis of the coil. 

Coil current 250 Amp. 

6. Computed lifting force and power adsorption .for a 
molten iron droplet along the axis of the coil. Coil 
current 300 Amp. 

7. Computed electromagnetic force field and computed 

6 3 

Joule heating (xlO W/ra ) in a levitated molten iron 
droplet. Coil current 250 Amp. 

8. Computed velocity field and temperature field 
in a levitated iron sphere. 

4 2 2 

9. Computed turbulent kinetic energy (xlO ,m /s ) RHS and the ■’ 

ratio of LHS. 

♦ 

10. Experimental carburization results (CO 2 .I 5 /CO 2 at 38.9 atm and 

165 0°c) (13,15) compared with the theoretical predictions 

for turbulent and laminar flow in the levitated drop. 

11. Computed velocity field and temperature field in a levitated 
molten drop assuming laminar flow. 


R13FKRENCES 


1. A.E, Jonkins, B. llairris and L.A, iiakcr, Symposium on 
metaliurcfy at hiqh pressures and hiqh temperatures, 

Met. Soc. AIME ConC. (1963), 22, p 23. 

2. E.C. Okress, D.M. Wroughton,’ G. Comenet55, P.H. Brace and 
J.e.R. kelly, J. Appl. Phy. (1952) , 23, 545. 

* 

3* J. Sisckelv and C.W. Chang, Ironmaking ans Steelmaking, (1977), 4, 

190. ‘ * 

4, J. Szokoly and C.W. Chang, Ironmaking and Steelmaking, (1977) 

4, 196. 

5. E. Tarapore and J. Evans, Met. Trans. B, (1976), 7B, 343. 

6. P, Cromer, Dr. Ing. Thesis, Universite de Grenoble, France, (1979). 

7, M. Choudhary and O'. Szekely, Met. Trans. B, (1900), IIB, 439. 

a: W.K.II. Panofsky and M. Phillips, Classical Electricity and 

Magnetism, 2nd edition, (1962), Addison-Wesley . 

9 . B.E. Launder and D.D. Spalding, Math. Models of Turbulence, 

Academic Press (1972) . 

10. B.E. Launder and D.B. Spalding, Computer Mcth. Appl. Mech. Bngr. 
(1974) , 3, 269. 

11. M. Choudhary, J. Szekely, B.I. Medovar, Y.G. Emelyanenko, 

Met. Trans. B. (1982), 13B, 35. 

12. N. E.-Kaddah and J. Szekely, j. Fluid Mech. 
in press. 

13. N.H. El-Kaddah, Ph.D. Thesis, Imperial Collpge, London, 

(1976) . 


ORIGINAL PAGE IS 
OF POOR QUALITY 


X<1. D.B. Spalding, Basic Equations of Fluid Mechanics, and 

Heat and Mass Transfer, and Procedures *for their Solution, 
Report HTS'17616, Imperial College, London, (1976). 

15. N.M. El~Kaddah and D.G.C, Robertson, Met.. Trans. E;, (197R) , 
9D, 191. 

16. D.G.C. Roberson and A.E. vTenkins, Heterogeneous Kinetics at 
Elevated Temperature, ed. Plenum Press, '(1.970), p 393 

17 . P. Bradshaw, Effects of Steamline Curvature on Turbulent 
Flow, (1973), AGAR Dograph 169. 

18. N.H.El-Kaddah and D.G.C. Robertson, Met. Trans. B, (1977), 

8B, 569. 

19. V. G. Levich, Physicochemicai Hydro“dynamics , Prentice-Hall 
Inc. (1962). 

20. J.T. Davies, Turbulence Phenomena, Academic Press. (1972). 

21. H. Branover and P. Gershon^ Symp. on Turbulent shear Flows, 
(1977) , Univ. Park Pennsylvania. 




ORIGINAL PAGE IS 
OF POOR QUALITY 


f 




% f 

1 







ORIGINAL PAGE 18 
OF POOR QUALITY 




Lifting Force FxlOf(N) 






omOINAL 

Of POOR QUALITY 



Distance from Center of Coil x lO^m) 


«— J.. ^ 


Power Adsorbed (W) 












Carbon (wt%) 





